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ABSTRACT 

The discovery of a population of high proper motion white dwarfs by Oppenhcimer et al (2001) has 
caused a lot of speculation as to the origin of these stars. I show that the age distribution of the white 
dwarfs offers a kind of sanity check in these discussions. In particular, this population appears to have 
a similar age distribution to those in the standard, thin disk white dwarf population. This is not what 
is expected for either the halo or thick disk, which are thought to be old populations. It may indicate 
a different dynamical origin entirely, or it may be just be the tail of a larger distribution, implying a 
significantly larger total density in even older white dwarfs. 

Subject headings: stellar dynamics — stars: kinematics — Galaxy: kinematics and dynamics 



1. INTRODUCTION 

The recent announcement of the detection of a signif- 
icant population of old white dwarfs with high proper 
motions (Oppenheimer et al 2001; hereafter OHDHS) has 
sparked a lot of interest in the nature of this population. 
If this population is indeed representative of the Galac- 
tic halo, as claimed by OHDHS, then it could represent a 
significant contribution to the claimed microlensing pop- 
ulation (Alcock et al 2000). However, the assignment of 
these white dwarfs to a particular dynamical family is still 
somewhat controversial. In particular, Reid, Sahu & Haw- 
ley (2001; hereafter RSH), by comparing with a local sam- 
ple of M dwarfs, suggest that this population is actually 
representative of the thick disk and consistent with known 
Galactic populations. In cither case, the OHDHS sample 
represents a new and interesting probe of the old stellar 
populations in the solar neighbourhood. 

Most of the recent controversy regarding the OHDHS 
claim centers on the dynamical interpretation of their re- 
sults; i.e. whether the detected white dwarfs represent 
the true halo or the thick disk and whether the density 
required is consistent with the expected stellar popula- 
tions or requires an interpretation in terms of a 'dark 
matter' contribution. In this article I wish to address an- 
other aspect of this population, namely the age distribu- 
tion and what that can tell us about the origin of these 
white dwarfs. First of all, in §|| I will review the argu- 
ments about the various dynamical populations suggested 
and the kinds of densities expected therefrom. Therafter, 
I will discuss in §|| the age distribution of the OHDHS 
sample and will discuss the implications in §|4| 

2. THICK DISKS & HALOS 

The OHDHS 'halo' sample of white dwarfs is culled from 
an area covering about 4000 square degrees towards the 
South Galactic Pole, or about ~ 10% of the sky. The 
magnitude limit of the search is R = 19.8, while the 
proper motion limit is 0.33".yr _1 . They apply a further 
velocity cut, selecting only those white dwarfs with ve- 
locity V± > 94km.s~ 1 (using photometric distance esti- 
mates) in order to define their high-velocity 'halo' popu- 



lation. Using a V max analysis they infer a mass density 
~ 1.3 x 10 _4 M Q .pc~ 3 in this high velocity population. 

This density is considerably larger than the estimate of 
1.3 x 10~ 5 M Q .pc~ 3 expected from the standard spheroid 
(Gould, Flynn & Bahcall 1998) and prompts interpreta- 
tions in support of halo white dwarf populations from mi- 
crolensing (e.g. Alcock et al 2000). However, alterna- 
tive explanations have already been suggested. RSH have 
shown that a similar velocity cut on their volume-complete 
sample of M dwarfs from Reid, Hawley & Gizis (1995) al- 
lows 20 out of 514 stars (~ 4%) into the 'halo' sample thus 
defined. As such, they interpret the OHDHS white dwarfs 
as simply a high-velocity tail of a disk population. This is 
also supported by the fact that the OHDHS 'halo' sample 
does show a distinctly asymmetric velocity distribution, 
suggesting some fraction belongs to a rotating component. 

To understand the various positions on this issue, we 
need to understand what is expected from each population. 
The thick disk is generally defined as a co-rotating (with 
the sun) population with vertical scale height ~ lkpc, ver- 
tical velocity dispersion ~ 40 — 45km.s~ 1 and a number 
density several percent of the local disk value (Sandage 
& Fouts 1987; Reid & Majewski 1994; Robin 1994). We 
should note, at this point, that the M-dwarf sample of 
Reid, Hawley & Gizis (1995) seems to contain a sig- 
nificantly larger fraction in the high velocity population 
(~ 15 - 20%). 

To infer absolute numbers, we review the range of Galac- 
tic mass models from Dehnen & Binney (1998). They as- 
sume a 5% fraction of the total disk mass in the thick com- 
ponent and find a range of total thick disk mass ranging 
from 8.3 x lO" 4 M .pc~ 3 to 1.11 x 10" 3 M o .pc~ 3 . How- 
ever, applying the velocity cut at 94km.s _1 removes all 
but the highest velocity objects from the sample. OHDHS 
claim this as a 2er cut. In reality, applying this cut to 
the proper motions drawn from a 3-D maxwellian only 
removes ~ 87% of the thick disk stars. Furthermore, al- 
lowing for a ~ 35km.s _1 lag in the thick disk rotation 
velocity (Chiba & Beers 2000), suggests as much as 20% 
of the stars make it into the sample. Thus, we may expect 
a range ~ 9.1 x lO- 5 M .pc- 3 -2.2 x lO" 4 M .pc" 3 . Hence 
the white dwarf density inferred by OHDHS appears to be 
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consistent with these numbers. However, one potentially 
disturbing feature of this agreement is that we have yet to 
account for the mass in stars below the turnoff i.e. those 
which have yet to form white dwarfs. 

To determine whether this is indeed disturbing, we have 
to infer the fraction of the total population mass that we 
expect to reside in white dwarfs. Assuming a Salpeter 
function from 0.1 Mq to 8M and assuming all stars with 
M > M TO = 0.82M Q form O.6M white dwarfs, the ratio 
of mass in white dwarfs to mass in stars with M < O.82M 
is 
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= 0.14 (M TO = O.82M ; x = 1.35) (3) 

so that ~ 13% of the total stellar mass is contained in 
white dwarfs. Thus, the expected total mass density in 
local, thick disk white dwarfs above the velocity cut is 
< 2.8 x 10~ 5 M Q .pc~ 3 , considerably less than the OHDHS 
determination (and similar to the Gould et al number for 
the spheroid). Of course, the Salpeter mass function di- 
verges at the low mass end, so this number could be sensi- 
tive to the low mass cutoff. Thus, let us repeat this calcula- 
tion using the empirical disk mass function of Gould, Bah- 
call & Flynn (1996) (and which doesn't diverge at the low 
mass end). Instead of a Salpeter slope x=1.35, this mass 
function has x=1.21 (note the different sign convention 
from Gould et al) between 0.6 - O.73M and x = -0.44 
for M < 0.6M Q . Gould et al argue that a correction for 
unresolved binaries will increase x to x ~ 0, so we adopt 
x — at the low end. Using this mass function (extrap- 
olated through the white dwarf region), the ratio of mass 
in white dwarfs to sub-turn-off stars is 0.41, i.e. a 29% 
share of the total mass budget. Thus, a more conservative 
estimate is 6.8 x 10~ 5 M Q .pc -3 thick disk white dwarfs in 
the OHDHS sample, i.e. only a factor of two smaller than 
the OHDHS number. 

Of course, there are several assumptions in such a naive 
model (mass function extrapolations, simple maxwellian 
velocity distributions) which can be changed to provide 
better agreement. We will consider these issues in ^J, but 
first we consider the other curious feature of the OHDHS 
sample, namely the age distribution. 

3. THE THICK DISK WHITE DWARF LUMINOSITY 
FUNCTION 

Again, it is useful to understand what we expect from a 
simple model. The thick disk is kinematically and chem- 
ically distinct (e.g. Freeman 1993; Majewski 1993) from 
the thin disk and is thought to be a population that formed 
primarily in a burst ~12 Gyr ago (e.g. Gilmore, Wyse & 
Jones 1995). Thus, we model this as a single burst of 
star formation 12 Gyr old, whose white dwarf ages vary 
depending on the mass and consequent main sequence life- 
time of the progenitor. To make a conservative model, we 
want to maximise the fraction of the white dwarf popu- 
lation that would be detectable by OHDHS. We assume 
a standard Salpeter IMF, a main sequence lifetime based 



on the models of Hurley, Pols & Tout (2000) and we shall 
assume that all stars above the turnoff mass (0.82 Mq in 
this model) make 0.5Mq white dwarfs (more massive white 
dwarfs cool more rapidly at late times because of earlier 
core crystallisation) , with hydrogen atmospheres (because 
these are the slowest cooling) and pure Carbon cores (be- 
cause these have the largest plausible heat capacity and 
thus cool most slowly). Thus, we are skewing the white 
dwarf luminosity function to the bright end. 

Figure [l] shows the resulting white dwarf luminosity 
function, using the models of Hansen (1999). For com- 
parison we include a luminosity function with all the same 
input parameters except that we assume a constant star 
formation rate over the last 12 Gyr. The most striking 
difference is that the burst population has a much sharper 
rise at the faint end. Thus, the vast majority (~ 90%) of 
the white dwarfs lie within 0.5 magnitude of the faintest 
white dwarfs in the burst case. This difference in scalings 
also demonstrates that one cannot simply try to scale the 
ratio of white dwarfs to M-dwarfs from the thin disk (a ~ 
constant star formation rate population) to the thick disk 
(if it is a burst population) as done in RSH. 

The next question to ask is whether the OHDHS white 
dwarfs are indeed old enough to account for the thick disk 
population. Again, we use our most conservative, 0.5 Mq, 
Carbon core, Hydrogen envelope models. Thus, we derive 
upper limits on the age. The absolute R-band magnitudes 
span the range Mr = 13 — 16. Using our conservative 
model, one finds an age of ~ 6.5 Gyr for Mr = 15 and 
~ 10.4 Gyr for Mr = 16. An age of 12 Gyr corresponds 
to Mr ~ 16.5, i.e. the oldest white dwarfs in the sam- 
ple are not older than 10 Gyr and most are considerably 
younger. Furthermore, recall that these are the most con- 
servative models, i.e. the ages are probably younger than 
this. Obviously the expected age of 12 Gyr is the combi- 
nation of both main sequence and white dwarf lifetimes, 
but, as we have shown above, we expect the white dwarfs 
to pile up at the faint end in a population resulting from 
a burst. Another way of putting this is to infer the main 
sequence mass of the progenitor by subtracting the in- 
ferred white dwarf age from the presumed 12 Gyr age of 
the burst. The faintest white dwarf (Mr ~ 15.7) has a 
white dwarf age < 9.9 Gyr (recall we are using the slowest 
cooling models) and thus comes from a star less massive 
than 1.35M if originating from a 12 Gyr burst. However, 
a more typical representative of the coolest OHDHS white 
dwarfs has Mr ~ 15.1, an age ~ 6.9 Gyr and a progen- 
itor mass ~ O.96M . If we adopt the Gould et al mass 
slope for the burst we find that (extrapolating to progeni- 
tor masses of 8Mq), we are missing the white dwarfs from 
between 50% and 80% (depending on the completeness of 
the OHDHS sample) of the total population. Since the 
observed densities are already slightly discrepant with our 
expectations, such enhancements of the density by factors 
of 2-5 would lead to a significant discrepancy. 

How much of this discrepancy could be due to inaccu- 
racies in the cooling models? Intercomparisons of differ- 
ent theoretical groups show model ages which differ by ~ 
10%, not enough to alleviate the above comparisons (a 
good review of the current state of the field can be found 
in Fontaine, Brassard & Bergeron (2001)). In particular, 
Fontaine et al present a model (pure Carbon core, Hy- 
drogen atmosphere) similar to that of our above conser- 
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vative model (although the mass is 0.6M Q ). We can com- 
pare the luminosity of such a 12 Gyr white dwarf from 
the two codes (and, of course, using a 0.6Mq model from 
our tables) and find logL/L Q = —4.73 (Fontaine) and 
logL/L Q = —4.78 (Hansen). The above white dwarfs are 
considerably brighter and younger. 

The most empirical measure of the ages of the OHDHS 
white dwarfs is to compare them directly to the thin disk 
sample of Liebert, Dahn & Monet (1988) ( although we will 
use the photometry of Bergeron, Ruiz & Leggett (1998)). 
Figure g shows that most of the OHDHS white dwarfs have 
similar ages to the thin disk white dwarfs, and a similar 
distribution. On the basis of this diagram alone, it would 
be difficult to distinguish this new population from the 
standard thin disk white dwarf population. 

Given that the OHDHS white dwarfs appear brighter 
than expected for their proposed parent population, one 
might wonder if there is some gross error in the photomet- 
ric distance indicator they used. The claimed 20% error 
is probably reasonable, given the small variation in radius 
between white dwarfs of different mass and composition. 
One potential note of caution though, is that mixed Hy- 
drogen and Helium atmospheres can wreak havoc with the 
colours of cool white dwarfs (e.g. the non-monotonic be- 
haviour of the mixed model shown in OHDHS Figure 4). 
This could contribute some additional systematic uncer- 
tainties. However, a wholesale overestimate of the dis- 
tances by this method is probably unlikely, given that 
many of the white dwarfs show Ha emission and the edge 
of that subset lies at an inferred Mr ~ 14.5 corresponding 
to temperatures ~ 5000 K (where we indeed expect the 
Ha to become undetectable). 

Finally, given these interesting results, one can ask what 
is required to see the kind of old white dwarfs we ex- 
pect. Searches in the R-band are optimal for the colours 
of old white dwarfs (Hansen 1998, 1999; Saumon & Jacob- 
sen 1999), and even the slowest cooling models described 
above predict Mr ~ 16.5 for a 12 Gyr white dwarf. To 
distinguish thick disk/halo white dwarfs we also require 
V_l > 100km. s -1 . Thus, we may set a target reduced 
proper motion H R = 23.1 + 5 log (V ± /100km.s _1 ) . The 
OHDHS sample is not complete at this level and thus there 
may still be a lot of white dwarfs to be discovered. 

4. CONCLUSIONS 

Based on the considerations of the previous section, it 
is tempting to conclude that the white dwarfs identified 
by OHDHS have a provenance similar to those in the thin 
disk. Certainly, on the basis of the colour-magnitude dis- 
tribution alone, it is difficult to distinguish the two pop- 
ulations. However, this sample was drawn from a high 



velocity sample. While arguments rage about their mem- 
bership in the halo or the thick disk, kinematic member- 
ship in the usual thin disk population seems unlikely. Of 
course, underlying this is the assumption that even old 
thin disk populations are described by maxwellian veloc- 
ity distributions with dispersions ~ 20km. s -1 . Perhaps 
this is simply indicating that there are other dynamical 
processes at work in the disk which can pump some stars 
up to velocities characteristic of the thick disk. 

The thick disk is generally regarded as kinematically 
and also largely chemically distinct from the old thin disk 
(Gilmore & Reid 1983; Sandage & Fouts 1987; Gilmore, 
Wyse & Kuijken 1989; Freeman 1993, Reid & Majewski 
1993). If the OHDHS white dwarfs are truly part of the 
thick disk they rather belie this distinction, suggesting 
that perhaps star formation in the thick disk is more com- 
plex than a simple burst at early times. One can make the 
OHDHS sample consistent with our preconceptions about 
the thick disk if one believes that it represents only the 
bright tail of a larger distribution. However, given the 
inferred densities as described in section a further ad- 
dition of significant mass would again imply a thick disk 
white dwarf contribution considerably larger than previ- 
ously thought. As such, it would imply a top-heavy mass 
function i.e. one more weighted towards the production 
of white dwarfs than a standard Salpeter mass function. 
Our naive estimates of §|| suggest as much as 80 % of the 
white dwarfs could lie beyond the OHDHS limits. Note, 
however, that the increased mass in white dwarfs would 
still not dominate the local disk mass, since the addition 
would only be of order the currently accepted (total) thick 
disk mass i.e. still an order of magnitude less than the 
thin disk contribution. 

Finally, much of the controversy regarding the interpre- 
tation of the OHDHS observations is motivated by whether 
this represents a detection of the population of stars re- 
sponsible for the microlensing observed by Alcock et al 
(2000). For this purpose it is not sufficient to distinguish 
between true halo and thick disk populations, as there are 
models (Gyuk & Gates 1999) which attempt to explain 
the microlensing results using rotating populations, remi- 
niscent of a very thick disk. Certainly, if one were to find 
significantly more white dwarfs (be they halo or thick disk) 
at even fainter magnitudes, the mass required would begin 
to point towards such a population. 

This work was supported by NASA through Hubble Fel- 
lowship grant HF-01120.01-99A from the Space Telescope 
Science Institute, which is operated by the Association of 
Universities for Research in Astronomy, Inc. under NASA 
contract NAS5-26555. 



REFERENCES 



Alcock, C. et al, 200, ApJ, 542, 281 

Bergeron, P., Ruiz, M.-T., & Leggett, S. K., 1997, ApJS, 108, 339 

Chiba, M. & Beers, T. O, 2000, AJ, 119, 2843 

Dchnen, W. & Binney, J., 1998, MNRAS, 294, 429 

Fontaine, C, Brassard, P. & Bergeron, P., 2001, PASP, 113, 409 

Freeman, K. O, 1993, in 'The Stellar Populations of Galaxies', IAU 

Symp. 149, ed. B.Barbuy & A.Renzini, p65 
Gilmore, G. & Reid, I. N., 1983, MNRAS, 202, 1025 
Gilmore, G., Wyse, R. & Jones, J. B., 1991, ApJ, 367, L55 
Gilmore, G., Wyse, R. & Kuijken, K., 1989, ARA&A, 27, 555 
Gould, A., Flynn, C. & Bahcall, J. N., 1998, ApJ, 503, 798 
Gyuk, G. & Gates, E., 1999, MNRAS, 304, 281 



Hansen, B. M. S., 1998, Nature, 394, 860 
Hansen, B. M. S., 1999, ApJ, 520, 680 

Hurley, J. R., Pols, O. R. & Tout, C. A., 2000, MNRAS, 315, 543 
Liebert, J., Dahn, C. C. & Monet, D. G., 1988, ApJ, 332, 891 
Liu, W. C. & Chaboyer, B., 2000, ApJ, 544, 818 
Oppenheimer, B., Hambly, N. C, Digby, A. P., Hodgkin, S. T. & 

Saumon, D., 2001, Science, 292, 698 
Reid, I. N., Hawley, S. L. & Gizis, J. E., 1995, AJ, 110, 1838 
Reid, I. N. & Majewski, S. , 1994, ApJ, 409, 635 
Reid, I. N., Sahu, K. & Hawley, S. L., 2001, /astro-ph/0104110 
Richer, H. B., Hansen, B. M. S., Limongi, M., Chieffi, A., Straniero, 

O., & Fahlman, G. G., 2000, ApJ, 529, 318 



4 



Robin, A. C, 1994, in 'Astronomy from Wide Field Imaging', IAU Sandage, A. & Fouts, G., 1987, AJ, 93, 74 

Symp. 161, ed. H.T.MacGillivray et al., p403 Saumon, D. & Jacobsen, S., 1999, ApJ, 511, L107 



5 




-5 I l l I I U I I 

o.oi io- 3 10- 4 10- 5 

L/L 



Fig. 1. — The lower solid curve is the distribution of white dwarf luminosities, assuming the model described in the text, for a single burst 
of star formation 12 Gyr ago. The upper dotted line is the distribution using the same model but assuming a constant star formation rate. 
The two distributions are normalised to have the same total number of stars formed. 
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Fig. 2. — We show here the absolute magnitudes (the 20% error bars on the photometric distances are not shown). The solid points are 
the OHDHS dwarfs (circles show Ha and triangles do not) and the open points are the thin disk dwarfs from Bergeron et al. The solid curve 
shows a O.5M0 pure Carbon core model with a Hydrogen atmosphere. Note that this is not intended as a fit to the observed sample but 
represents the slowest cooling white dwarf plausible. The positions of white dwarfs of ages 8, 10 and 12 Gyr are shown. 



